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FOREWORD 

The  major  goal  of  this  program  is  the  development  of  a breadboard 
in-house  laboratory  measurement  system.  In  pursuing  this  objective  many 
experimental  results  were  obtained  which  shed  light  on  some  improvements 
that  could  be  made  toward  further  development  of  such  a system.  These 
results  are  presented  in  this  report  along  with  some  tried  designs  and 
proposed  designs.  The  system  was  designed  to  readout  specklegram  data 
and  display  the  information  in  a readable  format  with  an  automated  opti- 
cal-digital set-up,  thus  reducing  the  analysis  time  greatly  over  existing 
methods.  The  overall  system  design  and  end  results  are  discussed  in  this 
report,  which  will  be  a final  report  to  the  Air  Force  Flight  Dynamics  La- 
boratory under  project  1367-02-12.  Mr.  Gene  Maddux  and  Dr.  Frank  Adams 
were  the  scientific  contacts  from  AFFDL. 

I would  like  to  thank  Dave  Flannery  of  AFAL/DHO-2  and  Chuck  Skinn 
of  System  Research  Lab.  for  their  assistance  in  pursuing  these  Program 
goals. 
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SECTION  I 
INTRODUCTION 

Near  the  beginning  of  FY74  Gene  Maddux  and  Dr.  Frank  Adams  of  the 
Air  Force  Flight  Dynamics  Laboratory  visited  Tom  Williamson  and 
Dick  Lane  of  the  Air  Force  Avionics  Laboratory  seeking  a feasible  method 
for  analyzing  specklegram  data  with  a semi -automated  optical  data 
processing  system.  A set-up,  made  up  of  existing  in-house  equipment 
including  a desk  calculator,  digital  voltmeter  A/D  converter,  and  a 
motor  driven  scanning  detector,  was  constructed  immediately  in  the 
Laboratory  of  Bldg  22B  and  after  several  test  conclusions  indicated 
that  an  optical  processing  system  could  be  fabricated  to  meet  the 
requirements  of  AFFDL.  (A  copy  of  the  initial  study  is  attached  at  the 
end  of  the  report.) 

In  January  1974  Dr.  Adams  requested  and  obtained  $24,000  of  AFFDL 
Laboratory  Director's  Funding  (LOF-74-1)  to  finance  a joint  program  with 
AFAL  for  developing  an  automated  data  readout  and  reduction  technique  for 
laser  speckle  photography.  AFFDL  funding  was  used  to  procure  equipment 
and  AFAL  provided  the  engineering  manpower. 

A plane  of  attack  was  then  begun  for  which  the  following  equipment 
was  ordered  at  the  beginning  of  FY75:  Minicomputer,  X-Y  Translation 
Stage,  2 dimensional  detector,  A/D  convertor,  clocking  circuit.  The 
final  delivery  of  equipment  was  received  in  May  of  1975  and  the 
fabrication  was  underway.  This  report  gives  the  demonstration  of  the 
brassboard  optical-digital  system  fabrication  for  performing  data 
reduction  of  specklegram  recordings. 
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SECTION  II 

STATEMENT  OF  PROBLEM 

In  the  early  1970s  several  researchers  in  the  United  Kingdom  showed 
that  small  in-plane  displacement  measurements  could  be  made  by  using  the 
inherent  speckle  properties  of  laser  light  and  employing  double  exposure 
photography  (References  1 to  6).  Further  development  of  this  technique, 
called  "Speckle  Photography,"  has  been  pursued  at  the  Air  Force  Flight 
Dynamics  Laboratory  (References  7 to  11).  It  was  found  that  the  speckle 
distribution  change  during  stressed  and  unstressed  conditions  on 
specified  surfaces  gives  a direct  indication  of  the  amount  of  surface 
motion  present  under  stress.  This  recording  process  is  performed  with 
a single  beam  laser  imaging  system,  and  as  many  as  40,000  data  points 
can  be  recorded  on  a 4 inch  by  5 inch  photographic  plate  with  one  double 
exposure.  At  this  point  it  was  realized  that  to  use  the  high  resolution 
recording  to  its  utmost,  a readout  system  had  to  be  established  to  speed 
up  the  method  for  extracting  the  data  in  its  entirety.  The  laboratory 
readout  method  used  originally  consisted  of  projecting  the  output 
fringes  on  a ground  glass  screen,  then  orienting  a slow  scanning 
detector  orthogonal  to  the  fringe  lines,  and  scanning  the  detector 
while  plotting  directly  onto  an  x-y  recorder.  This  process  has  to  be 
performed  ^'or  each  data  point  on  a specklegram  plate  with  the  addition 
of  one  accurate  method  for  scanning  from  data  point  to  data  point.  The 
average  point  to  point  recording  time  could  be  as  large  as  5 minutes, 
thus  requiring  140  hours  to  analyze  a single  plate. 

Careful  study  has  shown  that  to  establish  the  speckle  photography 
method  as  a feasible  tool  for  performing  non-destructive  testing,  the 
point  to  point  data  collection  time,  including  the  point  to  point  scan 
time,  must  be  1 second  or  faster.  At  these  rates,  the  maximum  time 
required  to  read  a specklegram  plate  at  maximum  resolution  would  be 
1.5  hours.  A faster  time  can  be  achieved  with  some  modifications  to 
the  existing  brassboard  described  in  this  report.  The  modifications 
are  discussed  in  the  section  entitled  "Conclusions  with  Alternate 
Solutions." 
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The  problem  of  data  reduction  of  specklegram  output  appears  to  be 
solvable  by  a digital-optical  system,  and  with  further  experimentation, 
a time  improvement  and  resolution  improvement  are  realizable  over  the 
initial  brassboard  setup. 
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SECTION  III 

DESCRIPTION  OF  OPTICAL  SETUP 

Figure  1 is  a graphical  illustration  of  the  optical  section  of  the 
brassboard  system  used  for  reading  out  specklegram  data  provided  by  the 
Flight  Dynamics  Laboratory  for  experimental  use.  A Spectra  Physics 
Model  124-A  He-Ne  Laser  with  15  milliwatts  of  power  and  a 1 millimeter 
diameter  beam  is  the  light  source  for  the  system. 

Two  lenses  with  a 4 to  1 focal  length  ratio  are  required  to  reduce 
the  laser  beam  diameter  to  its  optimal  operating  size.  Microscope 
objectives  of  20  power  and  5 power  were  chosen  for  this  experiment. 

The  lenses  are  oriented  such  that  the  distances  between  them  equals 
^20  * ^5’  ^ ™ beam  emerging  from  the  laser  enters 

the  specklegram  plate  at  .25  millimeters.  A .25  millimeter  beam  was 
chosen  because  experimentation  showed  that  a smaller  beam  would  create 
large  secondary  speckle  and  a largei  beam  would  reduce  the  inherent 
resolution  of  a data  plate. 

The  specklegram  plate  is  held  in  a plate  holder  attached  to  a 
5 inch  by  5 inch  Yoseniite  x-y  translation  stage  (Figure  2).  The  stage, 
with  .0001  inch  stepping  resolution,  mechanically  positions  the  data 
plate  to  each  new  data  point  position,  with  two  Slo-Syn  stepping  motors 
tied  directly  to  the  worm  gear  mechanism  on  the  stage.  Each  step  of 
rotation  by  the  motors  represents  a .001  inch  lateral  motion  of  the 
stage  and  is  controlled  by  the  PDP-11  central  processor  tied  to  a 
specially  built  electronic  circuit  constructed  in-house  by  AFAL.  The 
details  of  the  central  processor  software  and  control  circuit  are 
described  in  two  later  sections  of  this  report. 

A ground  glass  plate,  ground  with  10  micron  Carborundum,  serves  as 
a projection  screen  from  which  the  specklegram  fringes  are  imaged.  The 
ground  glass  screen  can  be  vibrated  as  a means  of  reducing  the  unwanted 
secondary  speckle  noise.  An  opaque  1/2  inch  diameter  spot  is  placed 
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Figure  2.  X-Y  Stepping  Stage 


on  the  ground  glass  at  the  center  of  the  speckle  pattern  canceling  the 
D.C.  component  of  the  readout  beam,  in  order  to  avoid  detector  saturation 
at  the  D.C.  imaging  point. 

The  primary  speckle  size  used  in  constructing  a specklegram 
determines  the  projection  cone  angle  of  the  fringe  pattern;  thus  as  a 
ground  glass  screen  is  moved  further  from  a specklegram  the  fringe 
size  increases  as  the  area  of  a plane,  subtending  a cone  perpendicular 
to  the  cone  axis.  When  sizing  the  secondary  speckle  to  match  the  ground 
glass  screen  consistency,  a hole  of  about  six  inches  in  diameter  is 
optimal;  thus  an  image  reducing  lens  is  required  to  project  the  proper 
fringe  size  onto  the  orthogonal  detector  assembly.  A vivitar  28 
Millimeter  focal  length  lens  is  used  for  the  brassboard  system  being 
described.  Since  specklegrams , as  stress  analysis  tools,  are  produced 
with  different  imaging  systems,  varying  reconstruction  cone  angles 
dictate  that  the  optical  geometry  between  the  specklegram  and  detector 
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be  variable.  The  ground  glass  screen  and  the  detector  assembly  are 
positioned  by  two  micrometer  controlled  slide  stages. 

Fulfilling  the  need  for  a medium  to  pick  up  the  imaged  fringe  data 
is  a custom  built  detector  system  consisting  of  two  linear  detector 
arrays  and  associated  operating  circuits  (Figure  3).  Appendix  V gives 
a detailed  description  of  the  specification  and  operating  procedures  of 
a Reticon  Corporation  Detector  System  similar  to  the  RL-64  WP  Model 
incorporated  in  the  optical  setup.  The  circuitry  for  operating  the 
detector  consists  of  five  main  electronic  sections:  1.)  - 2.)  Output 
drivers;  3.)  Clocking  circuit;  4.)  - 5.)  Sample  and  hold  circuits. 

The  output  drivers,  required  for  each  linear  array,  simply  amplify 
the  output  signal  from  each  detector  in  the  array  to  a 3 volt  maximum 
during  the  light  saturation  operation.  The  amount  of  light  required  to 
reach  saturation  is  obviously  dependent  on  the  length  of  time  the 
detector  is  receiving  light  and,  since  the  detectors  integrated 


Figure  3.  Reticon  Detector 
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all  the  incident  light  from  one  reading  to  the  next,  the  scan  time  from 
detector  to  detector  determines  the  optical  saturation  level.  The 
clocking  circuit  which  controls  both  detectors  has  a direct  influence 
on  the  detector  to  detector  read  time.  A more  detailed  description 
of  the  electronic  circuit  sequencing  will  be  discussed  in  Section  IV,  3. 
Circuits  4 and  5 are  merely  sample  and  hold  electronics  that  read  the 
output  from  the  driver  circuits  and  retain  them  for  a time  period  long 
enough  for  an  external  readout.  For  instance,  if  the  signal  illustrated 
in  Figure  4b  were  input  to  the  sample  and  hold  board  it  would  emerge 
as  shown  in  Figure  4a.  In  other  words,  a series  of  spiked  pulses  is 
transformed  into  a series  of  square  wave  pulses.  At  this  point 
attention  will  be  given  to  the  configuration  of  the  linear  detector 
element  relative  to  the  optical  fringe  output.  Figure  5 is  a magnified 
view  of  the  two  linear  detector  arrays  configured  in  a chip  holder.  The 
particular  geometry  chosen  will  allow  the  specklegram  fringes  to  be 
projected  onto  the  detector  array  at  any  angle  but  will  make  it  still 
possible  to  determine  the  orthogonal  spacing  between  fringes.  By  looking 
at  Figure  6,  we  see  that  if  the  x and  y axis  represent  the  two  linear 
detector  arrays  and  the  three  skewed  lines  represent  the  fringe  data, 
a simple  trigonometric  sequence  allows  us  to  determine  the  orthogonal 
spacing  (d)  and  the  fringe  angle  placement  (0)  from  the  fringe  spacing 
which  appears  on  the  detector  arrays  (x  and  y).  The  trigonometric 
sequence  is  as  follows: 

0 = tan  ^ y/x 
d = X sin  0 

d = X sin  (tan*^  y/x) 
d = xy/^y^+x^ 
d = X /s/(l  + (x2/y2)) 

Thus,  by  employing  two  linear  detector  arrays  set  in  a right  angle 
field,  a complete  set  of  useful  data  can  be  extracted  from  any 
specklegram  plate.  The  original  detectors  purchased  from  Reticon 
Corporation  were  64  elements  long,  and  since  a maximum  of  20  fringes 
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would  ever  enter  the  system,  the  Sampling  Theory  says  64  elements  would 
suffice.  But,  it  has  been  determined  through  experimentation  that  this 
number  is  too  low  to  perform  the  job  as  accurately  as  required. 


a 


b 


Figure  4.  Reticon  Output  Samples 


Figure  5.  Magnified  View  of  the  Reticon  Detectors  Elements 
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SECTION  IV 

DIGITAL  AND  ELECTRONIC  SETUP 

1.  COMPUTER  MODEL  PDP-11 -05 

The  nucleus  of  the  entire  specklegrarn  analysis  system  is  a Digital 
Equipment  Corporation  Model  PDP-11-05  minicomputer  with  4000  decimal 
word  storage  at  16  bits  per  word.  The  computer  performs,  under  program 
control,  all  calculations  on  the  incoming  data,  the  stepping  sequence 
of  the  x-y  translation  stage,  the  addressing  of  the  detector,  the 
scanning  of  the  detector,  and  the  retrieval  of  the  data  through  an 
A/D  converter.  Figure  7 is  a graphical  layout  of  the  processor  and  its 
peripheral  hookup.  A brief  sequence  of  extracting  data  from  a 
spe.klegram  might  be  as  follows:  1.)  Load  the  plate  into  its  holder 
and  position  it  to  a starting  point  at  the  edge  of  the  data  array. 

2. )  Set  up  the  optics  for  the  recorded  magnification.  3.)  Initialize 
all  peripheral  devices.  4.)  Enter  the  x-y  stage  point  differential 
(stepping  resolution).  5.)  Step  the  plate  to  its  first  data  point. 

6. )  Address  the  detector  array  and  perform  the  detector  scan. 

7. )  A/D  the  serial  data  into  core  memory.  8.)  Process  the  data  into 
fringe  spacing  and  angle.  9.)  List  the  data  to  the  teletype  and, 

10.)  Step  to  the  next  data  point  and  repeat  the  process.  At  the  end 
of  the  complete  process,  the  x-y  scan  parameters  can  be  re-entered  and 
a finer  or  coarser  scan  can  be  performed,  based  on  the  previous  data 
run.  At  this  point  a more  detailed  discussion  of  the  procedure  for 
performing  Steps  5 and  6 in  the  above  sequence  will  be  discussed. 

2.  X-Y  STEPPING  STAGE 

The  x-y  stepping  stage  (Step  5)  assembly  consists  of  a Yosemite 
5 inch  by  5 inch  translation  stage,  two  Slo-Syn  stepping  motors,  two 
Aerotech  Drive  Units,  and  a custom  built  plate  holder  mounted  to  the 
x-y  stage.  Stepping  of  the  motors  is  accomplished  by  shifting  pulses 
in  a serial  string  to  the  motors.  Each  pulse  rotates  the  motors  1.8° 
(200  pulses  - 360°),  and  the  polarity  of  the  pulse  determines  the 
direction  of  motor  rotation.  The  motor  driving  pulses  and  their 
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I 


ADC  (x) 


ADC  (y) 


Figure  7.  POP-ll  Computer  and  Peripheral  Hookup 
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polarization  are  generated  by  the  Aerotech  Drive  Units  and  only  signals 
inherent  to  or  input  to  these  drive  units  will  be  discussed  from  this 
point  on.  For  the  case  of  externally  driving  the  stage,  it  is  required 
to  apply  separate  5 volt,  50  visec  pulses  to  the  Aerotech  Drive  inputs 
for  rotating  the  motors,  and  either  a 0 or  5 volt  signal  to  the  direction 
inputs  to  dictate  positive  or  negative  direction  of  rotation.  Each  unit 
has  the  capability  of  running  internally  or  externally,  and  in  each  case 
the  number  of  pulses  can  be  monitored  by  a digital  readout  on  the  front 
panel  of  the  unit.  Figure  8 shows  the  2 Drive  Units  in  their  rack  mount. 
During  internal  operation,  for  use  in  initializing  plate  position,  the 
motors  can  be  stepped  in  single  pulse  mode  or  in  the  slow  mode;  this  will 
repetitively  transmit  pulses  to  the  motors  at  a rate  preset  by  a 
potentiometer  on  the  front  panel.  Another  necessary  feature  of  the 
internal  operation  is  the  ability  to  reset  the  digital  counter  at  any 
time  during  operation,  which  will  create  a zero  start  reference  point. 
Thus,  the  drive  units  can  be  set  to  internal  operation  for  initializing 
the  x-y  plate  position,  zeroing  the  digital  meter,  and  then  switching 
to  the  external  mode  for  the  minicomputer  to  take  control. 


I 


Figure  8.  Digital  and  Electronic  Hardware 
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Since  no  commercial  interface  cards  were  available  to  drive  the 
Aerotech  Units  externally,  a computer  interface  was  constructed  in-house 
to  generate  the  required  driving  signals.  Figure  9 is  a schematic  of 
the  interfacing  circuitry. 

Four  signals  had  to  be  generated  to  drive  the  x-y  stage  under 
computer  control,  two  directional  signals  and  two  scanning  signals. 

For  simplicity,  discussion  will  be  limited  to  only  the  x portion  of  the 
staging,  assuming  the  y portion  to  be  self-explanatory.  The  x-clock 
output  supplies  50  ;.sec,  5 volt  pulses  to  the  x-drive  unit  as  required, 
to  drive  the  x-scanning  motor.  In  order  to  get  a pulse  out  to  the 
x-clock  line  through  the  Figure  9 electrics,  a signal  has  to  be  applied 
to  lines  D1  and  El,  simultaneously  activating  the  D303  one  shot  gate, 
which,  depending  on  the  capacitor  value  tied  to  the  one-shot,  will 
produce  a low  voltage  50  usee  pulse  (capacitor  = .022  pf).  This  output 
pulse  is  then  amplified  to  a 5 volt  pulse  with  an  operational  amplifier 
tied  to  a 5 volt,  1 kilohm  resistive  netwoik.  The  problem  now  is  to 
determine  how  the  computer  uniquely  addresses  the  lines  D1  and  El. 

01  is  hooked  to  V2  and  an  output  from  this  AND  gate  is  achieved  by 
applying  a signal  on  BUS  DOO,  since  the  other  gate  input  is  always 
set  high.  The  three  BUS  lines  are  simply  binary  outputs  defined  as 
follows : 


LiiJ 

14  13  12 

11  10  9 

8 7 6 

5 4 3 

? 1 0 

f ^BUS  DOO 
•-BUS  DOl 
■■BUS  D02 


where:  the  numbers  1-15  represent  power  of  2, 

thus:  1 = BUS  DOO 

2 = BUS  DOl 

4 = BUS  D02 


El  is  hooked  to  the  AND  Gate  with  inputs  S2  and  Nl.  The  details  of 
these  outputs  can  be  referenced  in  the  Peripheral  Flandbook  under  the  M105 
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circuit  card  chapter.  For  all  practical  purposes,  in  understanding  this 
circuitry,  the  address  to  output  a signal  at  point  52  is  164000  and 
at  point  T2  is  164002.  In  total,  to  get  an  x-clock  output,  the  number  1 
has  to  be  moved  to  address  164000. 

In  a similar  manner,  a table  can  be  constructed  to  summarize  the 
four  x-y  stage  outputs  needed  for  computer  control; 


OUTPUT 

INPUT 

ADDRESS 

1 . 

x-clock 

o 

o 

II 

164000 

2. 

y-clock 

010  = 2 

164000 

3. 

x-di rection 

001  = 1 

164002 

4. 

y-direction 

010  = 2 

164002 

3.  RETICON  DETECTOR  ARRAY 

Since  the  specklegram  has  now  been  positioned  for  a data  point 
readout,  the  detector  has  to  be  addressed  to  read  the  fringe  data.  A 
brief  description  of  the  Reticon  Array  operation  would  be  of  assistance 
at  this  point.  For  simplicity,  the  discussion  will  be  limited  to  only 
one  of  the  linear  arrays,  because  the  electronic  inputs  affect  both 
arrays  simultaneously. 

Each  linear  array  consists  of  64  individually  addressed  detector 
elements,  which  are  interrogated  one  at  a time  by  a scanning  clock  and 
the  charge  accumulated  by  each  element  is  read  and  cleared  upon 
interrogation.  These  clock  pulses  can  be  generated  internally  by  a 
separate  clocking  circuit  at  a preset  rate  or  can  be  input  to  the 
clocking  circuit  from  an  external  source  at  any  rate  up  to  10  MHz, 
provided  the  pulses  are  at  least  .5  psec  in  duration  and  5 volts  in 
amplitude.  Whether  from  an  internal  or  external  clock  source,  the 
pulses  must  pass  through  additional  circuitry  to  be  counted  for  the 
purpose  of  generating  start  and  end  of  scan  pulses.  The  clocking 
circuit  board  has  12  binary  switches  (Figure  10)  for  setting  the  time 
delay  (number  of  pulses  delay)  desired  between  readings,  thus  the 
maximum  delay  is  2 -1  (8191  pulses).  At  the  end  of  each  scan  an  end 
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Figure  10.  Toggle  Switch  Configuration 


of  scan  pulse  is  generated  for  external  triggering  and  immediately 
before  each  scan  a trig  pulse  is  generated  also  for  internal  triggering. 
Referring  to  Figure  11,  a typical  scan  for  taking  a single  data  reading 
is  illustrated,  with  three  lines  of  information.  The  top  line  represents 
the  fact  that  the  data  is  being  read  to  the  computer  by  a negative  going 
pulse,  the  middle  line  shows  the  data  itself,  and  the  bottom  line  shows 
the  end  of  scan  pulses.  As  can  be  seen,  there  are  actually  two  complete 
scans  of  the  detector  elements,  before  the  third  scan  is  read  into  the 
computer.  The  reason  for  the  three  scans  is  that,  since  the  detectors 
integrate  the  light  incident  upon  them  up  until  they  are  interrogated, 
the  first  scan  most  likely  will  be  saturated  with  ambient  light  as 
seen  by  the  center  scan  at  the  far  left  of  the  figure.  The  second  scan 
sometimes  does  not  contain  all  64  elements  of  detector  information;  thus 
the  third  scan  is  the  first  complete  scan  containing  the  data  information 
only.  An  expanded  view  of  the  third  scan  computer  reading,  data  display, 
and  end  of  scan  pulse  is  given  in  Figure  12. 
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Figure  11. 


Figure 


Clock  Sequence  for  Read  Initialization 


. Clocking  Sequence  for  Read  Operation 
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Rpferrinc)  to  the  schematic  of  figure  9,  let's  now  look  at  the 
circuitry  built  to  generate  the  clock  pulse  and  A/D  the  data  to  the 
computer.  The  same  address  decoder  and  input  BUS  word  are  used  to 
generate  the  clock  pulses. 

OUTPUT  INPUT 

RETICON  100  = 4 

CLOCK 

A 200pf  capacitor  is  tied  to  the  the  D303  one  shot,  giving  a 5 psec 
pulse  output. 

Special  circuitry  had  to  be  constructed  to  match  the  Reticon  data 
output  to  the  input  of  the  4 channel  Analog  multiplexer  (Model  A-124). 
The  multiplexer,  A/D  converter  combination  requires  a 0 to  10  volt 

9 

input  swing  to  give  the  full  2 levels  of  digital  resolution;  thus  a 
variable  operational  amplifier  circuit  is  used  to  gain  the  maximum 
resolution  of  the  A/D  converter  by  taking  the  reticon  3 volt  signals 
to  the  10  volt  maximum.  The  end-of-scan  pulse  is  input  to  the  3rd 
channel  of  the  multiplexer  directly  and  is  used  to  synchronize  the 
software  program  to  the  detector  and  electronic  hardware. 

Two  cash  boards  had  to  be  attached  to  the  outputs  of  the  data  from 
the  reticon  to  act  as  a sample  and  hold  circuit  to  maintain  the  data 
from  each  detector  pulse  long  enough  for  the  A/D  converter  to  read 
its  value.  Figure  13  shows  two  data  outputs  before  and  after  the  cash 
board  circuitry.  Notice  the  short  length  of  time  the  detector  output 
remains  peaked  without  the  Cash  Board. 

Figure's  14,  15,  and  16  are  a few  examples  of  1,2,  and  9 fringe 
outputs  respectively.  In  the  lower  left  hand  corner  of  the  figure  is 
an  oscilloscope  trace  of  the  reticon  array  output,  with  its  analog  to 
digital  conversion  data  retrieved  from  computer  core,  to  the  immediate 
right  and  a graphical  plot  of  the  digital  data  is  displayed  above  the 
computer  printout.  Notice  how  closely  the  detector  output  matches  the 
graphical  plot  of  the  computer  data. 


ADDRESS 

164000 
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Figure  13.  Sample  and  Hold  Illustration 
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SECTION  V 
SOFTWARE 

1.  PDP-n  SOFTWARE  FUNCTIONS 

The  PDP-11-05  family  of  minicomputers  operates  on  simple  16  bit  machine 
language  with  7 registers  for  program  manipulation,  1 of  which  are  a stack 
pointer  (R6)  and  a program  counter  (R7).  The  stack  pointer  acts  as  storage 
pointer  for  keeping  track  of  subroutine  jumps  and  the  program  counter  keeps 
track  of  all  addressing  within  a program,  including  automatic  word  or  byte 
sequencing. 

The  method  chosen  for  writing  a program  depends  on  the  software  com- 
pilation routine  desired  by  the  user  to  generate  the  machine  language  for 
computer  operation.  The  highest  level  language  that  can  be  used  is  Fortran 
but  many  hardware  options  have  to  be  tied  to  the  computer  hardware  directly, 
and  along  with  the  high  cost  of  the  Fortran  compiler  software,  a substantial 
cost  increase  is  realized.  Due  to  the  limited  budget  available,  the  simplest 
of  all  languages  was  chosen,  and  other  than  creating  a substantially  more 
complex  programming  routine,  this  simple  language  proves  to  be  adequate  for 
doing  the  job.  This  simple  language  (assembly  language)  is  comprised  of 
the  following  programming  commands; 

1 . CLR  - clear 

2.  DEC  - decrement 

3.  INC  - increment 

4.  NEC  - negate 

5.  TST  - test 

6.  COM  - complement 

7.  ASR  - aritnmetic  shift  right 

8.  ASL  - arithmetic  shift  left 

9.  ADC  - add  carry 

10.  SBC  - subtract  carry 

11.  ROE  - rotate  left 

12.  ROR  - rotate  right 

13.  SWAB  - swap  bytes 

14.  MOV  - move 

15.  ADD  - add 

16.  SUB  - subtract 
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17. 

CMP  - 

comiiare 

18. 

RIS  - 

bit  set 

19. 

BIT  - 

bit  test 

20. 

BIC  - 

bit  clear 

21  . 

BR  - 

branch 

22. 

BEQ  - 

branch 

on 

equal 

23. 

BNE  - 

branch 

on 

not  equal 

24. 

BMI  - 

braiKh 

on 

minus 

25. 

BPL  - 

branch 

on 

plus 

26. 

BCS  - 

branch 

on 

carry  set 

27. 

BCC  - 

branch 

on 

overflow  set 

28. 

BVS  - 

branch 

on 

overflow  clear 

29. 

RVC  - 

branch 

on 

less  than 

30. 

BIT  - 

branch 

on 

greater  than 

31. 

BGE  - 

branch 

on 

greater  than  or  equal 

32. 

BLE  - 

branch 

on 

less  than  or  equal 

33. 

BGT  - 

branch 

on 

greater  than 

34. 

BHI  - 

branch 

on 

high 

35. 

RTS  - 

return 

from  subroutine 

36. 

JSR  - 

jump  tc 

1 subroutine 

37. 

HALT 

- halt 

38. 

WAIT 

- wait  for 

interrupt 

39. 

JMP  - 

j ump 

The  modes  available  for  applying  the  above  commands  are  listed  briefly  on 


the  following  tables  and  are  described  in  detail  in  the  PDP-11  Processor 
Handbook . 


1 MODE 

1 “ .1 

MODE 

NAME 

SYMBOLIC 

DESCRIPTION 

0 

REGISTER 

R 

R 

is  operand 

1. 

REGISTER  DEFERRED 

(R) 

R 

Is  address 

2. 

AUTO  INCREMENT 

(R)  + 

R 

is  address;  R + 2 R 

3. 

AUTO  INCREMENT  DEFERRED 

@ (R)  + 

R 

is  address  of  address;  R + 2 R 

4. 

AUTO  DECREMENT 

- (R) 

R 

- 2 ->■  R;  R is  address 

5. 

AUTO  DECREMENT  DEFERRED 

@ - (R) 

R 

- 2 -►  R;  R is  address  of  address 

6. 

INDEX 

X (R) 

R 

+ X is  address 

7. 

INDEX  DEFERRED 

0 X (R) 

R 

+ X is  address  of  address 

25 


* 


J 


AFAl  -TR-77-203 


MODE 

NAME 

SYMBOI.IC 

DESCRIPTION 

2 

IMMEDIATE 

n 

operand  n follows 

instr 

1 

Ali.MOI.nF 

(.)  A 

address  A follows 

instr 

h 

KK1.,\T[VE 

A 

instr  address  + 4 

+ X is 

address 

7 

REI..\TIVE  DEFERRED 

■1  A 

instr  address  + 4 

+ X is 

address 

of  address 

Ihri'o  examples  usinf>  the  commands  and  modes  listed  above  are 


^iven  as  follows. 
MOV  Rl.  (K21  + 


Ci.R  (R5) 


001 

000 

001 

010 

010 

mTiv 

MODE  R 

M(V)E  K 

010 

010 

111 

01 1 

100 

CMP 

MODE  7 

MODE  R 

LE 

000 

000 

100 

010 

1 10 

I) 

0 

0 

4 

2 

6 

Ll. 

000 

101 

000 

001 

101 

Cl.K 


MODE 
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2.  CONTROL  PROGRAM  SOFTWARE  DESCRIPTION 
I 


Due  to  the  limited  commands  and  the  lack  of  arithmetic  functions, 
the  software  required  to  control  the  complete  analysis  system  is  quite 
complex  and  lenqthy  contrary  to  what  miqht  be  expected  when  using  a 
higher  level  language.  The  program  can  be  outlined  into  21  different 
sections  or  subroutines  as  follows. 

1.  PROGRAM  START  AND  PARAMETER  LOADING  - START 

2.  ASC  II  AND  DECIMAL  TO  OCTAL  CONVERSION  - ASC  II 

3.  CONTROL  - MANIP 

4.  FAST  FOURIER  TRANSFORM  - FIT 

5.  DETECTOR  READ  - READ 

6.  DISPLACEMENT  & ANGLE  LIST  - PRINT 

7.  ANGLE  CALCULATION  - ANGLE 

8.  AMPLITUDE  CALCULATION  - AMP 

9.  TELETYPE  OUTPUT  - TTYOUT 

10.  OPERATION  SELECTOR  - OPER 

11.  TELETYPE  INPUT  - INPUT 

12.  DISPLACEMENT  CALCULATION  - DISP 

13.  FRINGE  EREQ  DETERMINER  - FREQ 

14.  OCTAL  TO  DECIMAL  CONVERSION  - OTD 

15.  X-Y  STAGE  CONTROL  - STAGE 

16.  PRINT  POSITION  - POSITN 

17.  DATA  BUEFERS  - RR0  - STEPLY 

18.  PRINT  BUFFERS  - ASKAV  - QUES 

19.  TRIG  TABLE  FOR  ANGLE  - TANBUF  - SINBUF 

20.  TRIG  TABLE  FOR  FFT  - COS  - SIN 

21.  SYMBOL  TABLE  - A - ZEROX 

Each  subprogram  can  be  described  briefly  in  the  order  that  it  appears 
in  the  overall  program  as  follows: 

1.  TITLE  merely  sets  a title  heading  for  the  program  that  appears 
at  the  top  of  each  page. 

5.  ENABLE  AMA  sets  the  assembler  to  define  the  machine  language  in 
the  Absolute  Memory  Addressing  mode. 
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8 - !i>  Those  lines  define  the  Idhels  assigned  to  the  bottofii  7 
words  of  core  in  computer  memory  with  register  numbers,  stack  pointer,  and 
program  counter. 

16  - 19  Lines  define  the  labels  assigned  to  the  listed  core 
locations  177560-177566  which  address  the  input  and  output  from  the 
teletype  terminal.  KBS  and  TPS  are  merely  status  words  and  KBB  and  TPB 
are  input  and  output  buffers. 

27.  500  is  the  starting  address  of  the  program. 

28.  START;  This  is  the  first  line  of  the  main  program  which  assigns 
a stack  pointer  address. 

30  - 31  Prints  "No.  of  Data  Averages"  by  addressing  the  TTYOUT 
Subroutine. 

32  - 37  Loads  ASC  II  value  of  No.  of  Averages  to  the  "MEAN"  buffer 
via  the  INPOT  subroutine  and  then  coverts  the  ASC  II  values  to  an  octal 
value  and  stores  it  in  the  TIMES  buffer. 

39  - 42  Clears  all  buffers  from  X LIMIT  to  STEP-X  in  the  data  buffer 
subroutine  called  RR0. 

43  - 75  Prints  enter,  and  convert  to  octul  : "X-Tot  Scan  Distance," 
"Y-Tot  Scan  Distance,"  "X-Increment,"  and  "Y-Increment,"  similar  to  the 
method  described  in  lines  32  - 37. 

68.  Sets  the  direction  of  the  X and  Y scan  to  be  positive.  Referring 
to  Figure  9,  move  a 4 to  164002  which  moves  zeros  to  164002  through  Bus 
DOO  and  DOl , making  the  output  to  the  direction  lines  0 volts  or  "low" 
produce  a positive  direction. 

76  - 77  Prints  the  header  title  as  shown  in  the  data  section  of 
this  report. 

78  - 79  Jumps  to  the  subroutine  "STAGE,"  which  initiates  the  X-Y 
stepping  stage. 

80  - 81  Prints  "X-Y  SCAN  COMPLETED." 

83.  Halts. 
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? 


86  - 9(1  ASC  11:  Suhtracts  60g  from  each  ASC  II  entry;  thus  the 
memher  6Ig(ASC  II  for  1)-6U^  = 1^^  and  71g{ASC  II  for  9)-60^  = Ilg,  which 
continues  until  a carriage  return  is  received  (<*I5). 

92  - 9A  Loads  Least  Significant  Digit,  which  is  already  in  octal 
format . 

98  - 104  Converts  2nd  digit  to  octal  format  as  follows: 

Add  12g  ^ duiiHiiy  register  for  each  unit  value  in  the  2nd  digit 

and  then  loads  the  final  value  to  R3,  which  is  the  finish  of  the  ASC  II 
routine,  and  is  the  final  octal  value  of  the  complete  list  of  digits. 

105  - 132  Converts  the  3rd,  4th,  and  5th  digit  in  the  same  manner 
as  in  steps  95  - 104. 

133  - 135  Stores  final  answer  in  OBUF. 

138  - 142  Clears  working  registers  and  buffers  from  RRO  to  KEEP 
(lines  892  - 911). 

143  - 140  Clears  FBUFAV  and  FFTBUF. 

147  - 153  Self  explanatory. 

154-158  If  the  frequency  is  below  4 fringes  then  prints  "LOW;" 
if  not,  then  prints  Go  To  a Frequency. 

159  - 162  Changes  Octal  Frequency  Count  to  Decimal  Value  and 
converts  to  ASC  II  Code;  then  print.svalue  in  Base  10  to  Teletype  Output 
Printer. 

163  - 175  The  first  time  through  MANIP  performed  complete  analysis 
on  X direction  of  Detector  data,  and  now  the  Y-direction  data  is  analyzed 
through  MANIP  and  printed  out  to  TTYOUT. 

178  - 462  FFT : Due  to  the  complexity  of  the  Fast  Fourier  Transform 
Routine  details  of  its  operation  can  be  obtained  from  Mr.  Lane  of  AFAL/DHO, 
and  will  not  be  covered  in  this  report. 

471.  READ:  Moves  a clock  pulse  from  the  Electronics  Unit  schematically 
illustrated  in  Figure  9 to  the  clock  input  of  the  Reticon  Array. 
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A72.  Moves  the  number  lOOSg  to  the  control  and  status  register 
(links  the  A/D  converter  to  the  software)  which  addresses  channel  2 of  the 
A/D  converter,  sets  the  priority  level  to  1 and  starts  the  A/D  operation. 

473.  Loads  the  CSR  (Control  and  Status  Register)  to  R0. 

474  - 475  Checks  for  the  Done  Bit  of  CSR  (Bit  7),  and  if  not  sets 
return  to  DONE  1,  otherwise  must  go  to  next  line  and  load  data. 

476  - 478  Moves  the  data  from  channel  2,  through  data  Register  176772, 
' to  register  RI.  Channel  two  is  hooked  to  the  detector  end  of  scan  (EOS) 

line  which  emits  a pulse  when  a complete  scan  of  the  64  detector  elements 
is  finished.  After  the  data  load,  the  program  checks  to  see  if  the  end 
of  scan  pulse  has  appeared.  If  not,  branches  to  CLOCK  1. 

479  - 482  Allows  for  2 EOS  pulses  before  continuing  to  the  next 
program  section. 

483  - 487  Continues  clocking  the  detector  through  the  dead  time 
between  scans  that  are  preprogrammed  through  the  toggle  switches  on  the 
clocking  circuit  board.  The  switches  are  set  to 

0000  nil  0011  = 227 

Thus,  it  works  out  that  202,  as  programmed  in  the  software,  provides  the 
proper  delay  to  allow  a data  read  at  the  correct  time. 
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488  - 481  Is  d software  delay  that  allows  a large  enough  time  between 
clock  pulses,  for  the  pulse  generator  of  Figure  9,  to  operate  properly. 

If.  for  instance,  the  MOV  = 4 , 0 » 164000  operation  were  performed  in  a 
small  program  loop,  the  pulses  generated  from  the  address  decoder,  arrive 
at  too  high  a rate  for  the  one-shot  delay,  thus  exiting  at  the  clock 
pulse  output  as  a voltage  level,  not  a pulse. 

492  - 494  Decides  whether  to  read  the  X detectors  or  the  Y detectors. 

495  - 500  The  entire  clocking  routine  previous  to  these  steps  were 
initialization  procedures  aimed  at  setting  up  the  A/D  converter  to  read 
data  from  the  detector  array  at  a time  when  meaningful  data  was  present. 
These  few  lines  read  this  unambiguous  data  and  load  it  into  core  memory. 
{R4)  is  the  storage  medium  where  R4  points  to  either  FFTBUF  or  FBUFAV. 

505  - 512  If  averaging  of  the  data  is  desired,  which  is  a buffer 
that  designates  a data  load  to  either  FFTBUF,  or  FBUFAV. 

513  - 525  Repeats  the  data  read  operation  for  the  Y-axis  with  the 
same  logic  as  in  steps  492-509. 

526  - 540  Averages  the  data  and  loads  the  averaged  data  to  the 
FFTBUF  buffer.  Since  the  averaging  method  probably  does  not  appear 
obvious,  a brief  description  follows:  523  - 525  TIMES  contains  the 
number  of  desired  averages  with  TRACK  incrementing  with  each  average 
performed  (TRACK  = TIMES  means  done).  526  - 527  The  data  is  loaded  to 
FFTBUF  from  the  first  data  read;  528  - 534  after  first  two  data  reads, 
each  data  buffer  (FFTBUF:  FBUFAV)  is  divided  by  2 (ASR)  added  together 
and  loaded  to  FFTBUF.  535  - 537  Prepares  to  load  new  data  to  FBUFAV 
for  another  average. 

543  - 549  Allows  for  reduction  of  redundancy  between  first  and 
second  half  of  the  detector  active  elements  by  loading  the  data  stored 
for  the  first  half  of  the  detector  into  the  storage  location  delegated 
for  the  second  half. 

555.  PRINTS. 

556  - 562  Clears  the  XDISTO  and  YDISTO  buffers;  then  calculates  the 
X-spacing  from  the  X-frequency  as  follows:  (f/No.  of  detectors  = spacing) 
where  f " ROBUF,  No.  of  detectors  = 64.|q  = 100^^. 
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563  - 567  Same  as  above  for  Y-spacing. 

568  - 576  Goes  to  Angle  calculation  subroutine;  then  returns  and 
prints  answer.  Goes  to  Displacement  subroutine;  then  returns  and  prints 
answer. 

581.  ANGLE. 

582  - 593  Loads  X and  Y spacinqs;  multiplies  the  Y spacing  by  lOOg 
to  obtain  a whole-number  output  from  the  division  of  Y by  X.  The  angle 
is  calculated  by: 

0 = tan  Y/X 

594  - 599  Takes  the  value  Y/X  (stored  in  RO)  and  scans  it  through 
the  tangent  table  (TANBUF),  listed  in  one-degree  increments,  until  the  two 
values  compare;  then  stores  the  comparison  address  to  RATIO. 

600  - 607  Clears  DBUF  then  Jumps  to  Subroutine  OTD  and  loads  to  DBUF. 

612.  AMP. 

613  - 620  Loads  start  address  for  Real  Part  of  FFT  (FFTBUF)  and 
imaginary  Part  of  FFT  (FFTBUF  + 200).  Checks  to  see  if  data  point  of  real 
and  imaginary  are  negative,  if  so  make  positive. 

621  - 628  Squares  the  real  data  point  amplitude. 

629  - 641  Squares  the  imaginary  data  point  amp. 

637  - 641  Adds  real  and  imaginary  data  points  and  moves  to  original 
real  data  point  storage  location  and  checks  to  see  if  all  data  points  have 
been  squared;  if  not,  increment  the  FFTBUF  addressed  for  real  and  imaginary 
and  repeats  the  complete  process  starting  at  line  615. 

651.  TTYOUT. 

651  - 656  Tests  the  Teleprinter  Status;  if  read,  moves  data  to 
Teleprinter  Buffer  and  prints.  Checks  if  finished  printing  by  a 0 data 
point  appearing. 

661.  OPER. 

661  - 678  Checks  if  a given  letter  (S,  X,  Y,  or  0) 
was  input  from  the  Teletype  and,  if  so,  load  R5  accordingly. 
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683.  INPUT. 

683  - 687  Checks  the  Keyboard  Status  and,  if  ready,  moves  the  ASC  II 
code  of  the  desired  key  to  R4;  then  clears  the  part  bit  (Bit  7)  and  moves 
the  answer  to  address  of  (RO). 

696  - 702  Prints  correct  character;  checks  if  input  is  completed 
{=15  - Carriage  Return)  and  moves  a line  feed  to  the  printer  (*12). 

704.  DISP. 

704  - 707  Moves  Sin  Buffer  address  to  R1  and  adds  RATIO  to  it. 

Ratio  contains  the  address  representation  of  the  angle  calculated  in  the 
ANGLE  subroutine.  Since  the  Displacement  is  calculated  by 

d = X sin  0 

and  the  sin  of  '■  can  be  calculated  from  a table  similar  to  TANBUF,  called 
SINBUF,  which  is  also  set  up  in  one-degree  increments;  d is  easily 
obtained  by  multiplying  X by  the  value  obtained  from  the  SINBUF  table. 

708  - 723  Multiplies  X DISTO  by  sin  0 and  divides  the  answer  by  lOOg 
since  SINBUF  value  is  premultiplied  bv  lOOg;  then  the  answer  is  changed 
to  decimal  and  control  is  returned  to  PRINT. 

727.  FREQ. 

729  - 731  Checks  for  a high  amplitude  D.C.  component  in  the  FFT 
Data;  if  lower  than  156,  sets  HIDC. 

732.  FFTBUF  + 140  represents  a fifteen  Hertz  frequency  bin  and, 
for  determining  the  frequency,  it  is  assumed  that  no  frequency  will  be 
higher  than  15  Hz  or  15  fringes. 

733  - 738  Checks  for  highest  amplitude  frequency  bin  by  comparing 
the  amplitude  of  the  presently  addressed  bin  to  the  maximum  amplitude  of 
all  the  previously  searched  bins;  if  larger  goes  to  HOLD.  If  it  is  smaller, 
increments  the  next  bin  and  either  returns  to  search  at  SRCH  or  moves  to 
QUIT  if  all  bins  have  been  searched. 

739  - 747  For  HIDC,  checks  data  with  #10  amplitude;  if  below 
GTO  QUIT  - 6. 
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748  - 7S0  If  froquency  was  valid,  loads  address  to  R2  and  the 
amplitude  to  R3,  the  increment  between  data  point  scans. 

753  - 757  If  no  frequency  was  stored  (R3  - 0),  then  moves  1 to  NONE. 

758  - 765  Looks  at  amplitudes  on  both  sides  of  the  highest  amplitude 
fretiuency  bin  (R3)  and  determines  if  the  actual  frequency  is  slightly 
higher  or  lower  than  R3  by  loading  1 to  R5,  if  above;  and  loading  0 to  R5, 
if  below. 

766  - 769  Since  FFTBUF  + 200  is  0 frequency  and  each  bin  in  the 
FFTBUF  is  one  cycle,  then  subtracting  R2  from  FFTBUF  + 200  and  dividing 
the  answer  by  two,  results  in  the  final  frequency  in  Hertz. 

770  - 777  Loads  the  frequency  to  OBUF  and  a + or  - to  OBUF  + 2. 
Returns  control  MANIP. 

786.  OTD. 

786  - 822  Changes  data  in  OBUF  to  a base  10  equivalent;  then 
converts  to  ASC  II  code  for  printout. 

827.  STAGE. 

832  - 840  MOV  #1,0#  164000  sends  a X-clock  pulse  to  the  X-Y  stage 
driver.  Unit  keeps  track  of  the  number  of  pulses  transmitted  for  each 
scan  between  data  points,  whereas  LINE  keeps  track  of  the  total  number  of 
pulses  transmitted  to  the  total  X-scan  for  one  line.  DELAY  simply  stores 
a programmed  delay  to  allow  the  pulses  to  zero  before  another  signal 
appears.  Unit  is  compared  to  STEP  X which  contains  the  data  determining 
the  increment  between  data  point  scans. 

845  - 848  Checks  to  see  if  the  end  of  line  distance  is  reached; 
if  so.  Branch  to  2ER0X. 

865  - 873  Changes  the  sign  of  X and  sends  pulses  to  drive  the  X 
unit  until  the  X position  is  0 or  its  original  position.  A long  delay 
is  required  (40000)  before  further  execution  of  scanning  to  allow  for 
experimental  stabilizing.  At  this  point  the  X sign  is  changed  again  and 
a Y scan  is  performed. 
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883.  805 1 TN. 

I 

883  - 886  Clears  the  Decimal  Buffer  (DBUF). 

887  - 890  Loads  the  X-position  in  octal  and  converts  to  decimal 
for  terminal  printout. 

891  - 896  Same  as  above  for  Y-position. 

897  - 900  Reloads  the  data  point  scan  separation  and  returns  to 
STAGE  subroutine. 

The  last  two  pages  of  the  program  are  a table  showing  the  numerical 
address  assigned  to  all  Labels  in  the  program. 

The  description  of  the  software  for  operating  the  overall  System, 
has  been  explained  only  in  brief  for  the  purpose  of  this  report  and  a 
more  detailed  discussion  can  be  planned  upon  request  by  the  Flight 
Dynamics  Laboratory. 


t 
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SFCTION  VI 
DATA 

In  this  section  of  the  report  we  will  discuss  the  individual  data  shown 
in  Tiqure  17  throuqh  30  and  the  siqni f icance  of  the  overall  data  collection 
as  presented  in  the  Results  Section.  A complete  data  collection  operation 
had  not  been  perfoniied  at  the  writinq  of  this  report,  but  shortly  there- 
after, it  was  anticipated  that  a complete  set  of  data  would  become 
available  for  listinq  to  the  Fliqht  Dynamics  Laboratory. 

The  first  method  pursued  for  analyzinq  speckleqram  data  was  to  look 
for  amplitude  peaks  or  nulls  in  the  detector  scanned  data  siqnals  and 
averaqe  the  distance  between  each  peak  or  null  by  countinq  the  number  of 
detectors  lyinq  between  them. 

for  initial  experimental  trials  a set  of  fringe  neqatives  were 
constructed  photoqraphical ly  to  input  to  the  optical  systems  as  a sample 
of  high  contrast  fringes  for  setting  up  the  computer  software.  Some  data 
from  these  initial  runs  is  shown  in  Figures  17,  18,  19,  and  20.  Frequencies 
of  5 Hz,  6 Hz,  7 Hz,  8 Hz , 9 Hz,  10  Hz,  20  Hz,  30  Hz,  and  40  Hz  were  the 
photogranhic  inputs.  As  can  be  seen,  the  system  operated  nearly  flawlessly 
for  these  types  of  fringe  inputs.  At  this  point  some  real  fringe  data  was 
input  to  the  optical  system  and  a couple  of  the  best  data  runs  are  shown 
in  Figures  21  and  22.  Notice  the  variance  in  what  appears  to  be  peaks  or 
noils  in  the  digital  data  printout.  The  data  for  most  cases  was  less 
impressive  than  that  shown,  thus  a new  method  for  counting  the  fringe  data 
had  to  be  establ ished. 

Coincidental  with  our  search  for  a new  method  of  analysis,  an  electronics 
equipment  oriented  company  made  available  an  FFT  (Fast  Fourier  Transform) 
operator  for  demonstration.  This  turned  out  to  be  the  opportunity  we  needed 
to  see  if  an  FFT  program  would  solve  some  of  the  problems  we  had  encountered 
during  peak  detection.  A typical  speckleqram  data  was  projected  onto  a 
Reticon  Array  that  we  supplied  and  input  to  the  FFT  hardware  they  supplied. 
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Figure  17.  Data  Runs  With  High  Contrast, 
Ideal  Fringe  Pattern 
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Figure  18.  Data  Runs  With  High  Contrast, 
Ideal  Fringe  Pattern 
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Figure  20.  Data  Runs  With  High  Contrast 
Ideal  Fringe  Pattern 


Figure  19.  Data  Runs  With  High  Contrast, 
Ideal  Fringe  Pattern 
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Figure  21.  Data  Runs  With  Real  Specklegram  Data 
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Figure  22.  Data  Runs  With  Real  Specklegram  Data 
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Fi<)tire  23  is  ci  four  plot  sequenced  illustration  of  the  result  of  this  mock 
setup.  The  first  plot  shows  the  specklegram  data  after  and  before  averaging 
512  times  (top  and  bottom  respectively).  The  second  plot  shows  the  first 
TFT  with  a strong  peak  to  the  far  left  and  a secondary  peak  to  its  injiiediate 
right.  The  first  peak  turned  out  to  be  the  reticon  clock  frequency  and  the 
second  peak  the  actual  data  frequency.  The  third  plot  shows  the  clock 
frequency  filtered  out,  thus  rendering  the  strongest  frequency  peak  to  be 
4.2  KHz  instead  of  the  200  Hz  Reticon  clock  frequency.  Plot  four  shows 
the  inverse  transform  of  the  FFT  data  giving  back  a signal  with  the  same 
frequency  as  plot  number  1.  Thus  it  was  concluded  that  the  FFT  approach 
would  become  the  solution  to  the  problem. 

After  an  extensive  amount  of  prograimiing  was  performed  to  implement 
an  FFT  software  routine,  the  photographic  fringes  were  introduced  into  the 
system.  Figures  24,  25,  26,  27,  and  28  show  a data  run  with  different 
frequencies  introduced  into  the  x and  y detectors  of  the  Reticon  Array. 
Notice  please,  that  each  frequency  calculated  is  as  nearly  right  as  can  be 
interpreted  from  the  Reticon  detector  plot  shown  on  the  scope  traces.  The 
pluses  and  minuses  listed  before  each  frequency  indicated  whether  the 
frequency  displayed  is  slightly  higher  or  lower  than  the  frequency  printed. 

The  computer  listing  of  Figure  29  is  a typical  teletype  printout 
obtained  during  a data  run  with  the  software  advanced  as  far  as  it  is 
today.  The  first  question  asks  if  the  user  wants  a list  of  the  FFT  data. 

The  second  question  asks  how  many  times  he  wants  the  fringe  data  input 
and  averaged  before  performing  the  FFT.  "OUTPUT"  asks  if  an  x or  a y 
detector  scan  is  desired.  The  next  few  questions  set  up  the  scanning 
parameters  for  the  x-y  scanning  stage,  and  are  self-explanatory.  At 
this  point  the  computer  takes  over  and  prints  out  the  x and  y position  of 
the  data  plate  in  the  x-y  translation  stage,  the  calculated  value  for  the 
frequency  of  the  fringes,  and  a "LIST,"  which  is  merely  a listing  of  the 
Fourier  Transform  of  the  averaged  data,  so  as  to  determine  if  the  cal- 
culated frequency  agrees  with  the  Transform  data. 


I 
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Figure  23.  Sample  Data  Analyzed  With  Hardware  FFT 
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OUTPUT  ? :: 

F '.SQUaiCVs  -830  10 


OUTPUT  ? V 
F r.i:,GUU;JCV=  +83003 


OUTPUT  ? 


Figure  24.  Data  With  Software  FFT 


OUTPUT  ? X 
F :.£UU£!iCV=  +30005 


OUTPUT  ? Y 
FKEeU£NCY=  +00008 


OUTPUT  ? 


Figure  25.  Data  With  Software  FFT 
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OITP'UT  ? X 

f r.'£QUh;r;cv=  ♦e00f^e 


OVTPUT  ? YY 
f r.£.'.lUi.tJCY=  -B0K08 


Figure  26.  Data  With  Software  FFT 


OUTP’JT  ? X 
Fr.iwU^jCY=  +00006 


OUTPUT  ? Y 
f .PEOUZ.NCY= 


• 0 300  6 


OUTPUT  ? 


Figure  27.  Data  With  Software  FFT 
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LISTED  ? Y 

NO.#  DATA  AVETJ^GES  ^ 1 

OUTPUT  ? Y 

X-TOT  SCAi'I  DISTANCE  ? / 1000 

Y-TOT  SCAN  DISTANCE  ? , 300 

X-INCr.EMENT  ? , 200 

Y-IHCP.EMENT  ? , 100 


X Y FREQ  LIST 

00200  00000  LOU  226>00U00I,001,035,001>002>001, 

00 1, 00 U 00  1 >000, 034^030, 00  1, 

00i:i00  00000  0002  122,025,002,300,002,031,002,003, 

001.300.001.003.304.000. 031, 

00630  00000  0007  065,009,010,330,301,001,026,039, 

021.000. 002.001.331.001.000, 

00803  00000  0008  065,001,002,001,301,030,002,009, 

001.001.032.000. 031.930.090, 

0 1000  00000  j;007 _0  50,  00  1 , 00  1 , 204,  00  1 , 03  1 , 00  5,  034, 

004,000,001,001,001,033,004, 

00200  00100  LOU  200,094,331,004,091,901,022,029, 

000, 00 1,0 3 1,390, 0 3 5, 3 0 9, 0 00, 

00200  00100  LOU  197,000,001,901,302,030,931,001, 

004.099.002.001.004.090.031, 

00400  00100  0002  170, 025, 00 1 ,999,  902, 20 ;, 202, 920, 

‘ 001, £01, 005, 003,301, 000, 301, 

00600  00100  0007  050,901,991,094,001,039,337,923, 

005.030.990.001.001.900.031, 

0 0800  00100  000  5 0 65, 991, 00  5, 001, 01  7,  0 91, 005, ':’91, 

00  5,030,0.02,001,091,091,031, 

0 1090  00100  0006 0 50,  094,  002,  30  1 , 30  2,  099, 305,  3.P4, 

005,001,0£5,090,905,0GO,09t, 

00000  90200  LOU  122,002,002,091,002,991,002,091, 

004,002,001,300,004,009,091, 

90200  00200  LOU  170,005,001,001,002,931,309,391, 

005,002,994,001,001,030,001, 

00400  03200  LOU  _ 170,004,301,330,892,031,C32,C0l,  .. 

00 1.0 09. 092. 000.  090, £30, 031, 

00600  00200  0005  085,004, 091, 001, 917, 016, 332, 901, 

002.000. 005.090.091.331.339, 

00800  00200  0006  032,00 1 , 00  5, 033, 03  5, 935, 00  I , 030, 

005.001.035.000. 009.801.331, 

01000  00200  0007  082,901,802,004,891,900,039,001, 

004,001,880,004,804,0  30,90  1, 

00900  00300  0011  082,094,035,098,081,091,001,331, 

001.000. 805.008.001.001.000, 


Figure  29.  Finalized  Data  Run.  Sample 
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Figure  30;  a-r  gives  a comparison  between  the  data  incident  on  the 
detector  and  the  data  plots  from  the  exit  lines  of  the  scanning  detector, 
monitored  on  an  oscilloscope.  The  top  line  of  each  plot  is  a detector 
scan  in  the  x-direction  and  the  bottom  plot  is  a detector  scan  in  the 
y-direction.  It  is  left  to  the  reader  to  interpret  the  results. 


47 


AFAL-TR-77-203 


I 


I 


I 


t 


I 


Figure  30  (a-c).  Real  Data  With  Different  Specklegram  Outputs 
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Figure  30  (g-i).  Real  Data  With  Different  Specklegram  0..tputs 
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Figure  30  (j-1).  Real  Data  With  Different  Specklegram  Outputs 
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Figure  30  (m-o).  Real  Data  With  Different  Specklegram  Outputs 
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Figure  30  (p-r).  Real  Data  With  Different  Specklegram  Outputs 
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SECTION  VII 
RESULTS 

The  results,  although  they  did  not  meet  the  original  regui rements , 
produced  some  meaningful  information  as  well  as  some  aid  toward  develop- 
ment of  a system  to  solve  some  of  the  problems  of  specklegram  readout  for 
the  Air  Force  Flight  Dynamics  Laboratory. 

The  discussion  that  follows  in  this  section  pertains  only  to  the 
system  as  it  exists  in  its  present  state  of  development,  directing  itself 
in  a time  ordered  fashion  from  the  program  beginning  to  now. 

Much  time  was  spent  at  the  start  in  generating  the  correct  electronics 
control  signals,  through  the  specially  built  logic  circuitry,  for  activating 
the  x-y  stepping  stage  and  the  two  dimensional  detector  array.  Many 
capacitor  and  resistor  changes  were  experimentally  adjusted  to  allow  the 
proper  amplitude  and  pulse  widths  of  the  driving  pulses  required  to 
operate  these  peripheral  units  optimumly. 

Almost  fifty  hours  of  machine  shop  time  was  required  to  build  the 
hardware  to  integrate  much  of  the  equipment.  For  instance,  a specklegram 
holder  had  to  be  built  for  mounting  to  the  x-y  stage,  a base  was  constructed 
to  mount  the  x-y  stage  to  the  optical  bench,  a base  was  bu'lt  for  mounting 
the  detector,  and  some  special  hardware  was  constructed  for  rack  mounting 
the  electrical  circuitry  and  labeling  the  faceplate. 

The  system  was  then  assembled  and  checked  out  in  its  final  form  as 
satisfactory  and  operable. 

The  software  required  most  of  the  time  from  this  point  on.  Each 
subroutine  discussed  in  the  software  section  of  this  report  was  written 
and  rewritten  until  a final  workable  program  was  achieved.  The  data 
analysis  section  of  the  software  program  was  written  in  two  different 
forms;  1.)  peak  detection,  2.)  Fast  Fourier  Transform.  Essentially 
all  the  software  for  the  peak  detection  was  thrown  out  after  it  was 
determined  to  be  inadequate  for  data  analysis.  Nearly  all  of  the  software 
progressions  for  both  methods  will  be  available  at  the  contract  close-out 
time. 
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SECTION  VIII 
CONCLUSIONS 

In  this  section,  I would  like  to  discuss  the  present  capability  of 
the  final  system  and  also  interject  some  possible  solutions  toward 
establishing  a more  accurate  and  feasible  system  for  carrying  out  the 
initial  objectives  laid  out  for  this  project. 

It  appears  that  most  of  the  objectives  such  as  data  analysis  speed, 
automatic  data  input  with  .001-inch  accuracy,  and  fringe  angle  and 
spacing  detennined  by  nonmechanical  means,  have  been  achieved  with 
relatively  few  problems.  It  has  been  demonstrated  that:  1.)  a data 
analysis  rate  of  1-2  Hz  is  feasible  in  the  existing  software  system; 

2.)  a .001-inch  accuracy  of  stepping  positions  can  be  achieved  under 
computer  control  at  the  above  data  analysis  rate;  and  3.)  the  detector 
shceme  of  two  orthogonal  linear  arrays  allows  for  the  determination  of 
angle  and  displacement  completely  independent  of  the  fringe  orientation 
on  the  detectors. 

After  completion  of  the  software  preliminaries  and  a checkout  of  the 
FFT  portion  of  the  computer  program,  some  of  the  specklegram  data  furnished 
by  Flight  Dynamics  Laboratory,  was  introduced  for  initial  tests.  During 
the  testing  period,  various  optical  configuration  changes  were  made  such 
as  varying  the  size  of  the  readout  laser  beam  and  changing  the  image  size 
of  the  specklegram  fringes  on  the  detector  faces,  which  compensate  for 
noise  and  detector  optimization  respectively.  It  was  found  that  a beam 
size  of  approximately  .25  mm  was  optimum  to  give  maximum  fringe  contrast 
yet  not  cause  too  high  a noise  level  due  to  large  secondary  speckle.  It 
was  also  determined  that  an  image  size  about  1-1/4  times  the  detector 
length  utilized  the  detectors  to  their  fullest,  maintaining  the  best  ratio 
between  number  of  detectors  per  fringe  and  number  of  fringes  per  detector. 
For  example,  if  15  fringes  were  present  across  the  complete  halo,  then 
approximately  12  fringes  would  be  imaged  across  the  64  detector  elements; 
thus  there  would  be  5-1/3  detectors  per  fringe.  For  a low  number  of 
fringes  (5)  there  would  be  4 fringes  projected  onto  the  detectors  thus 
16  detectors  per  fringe. 


55 


AFAL-TR-77-203 


The  best  results  that  could  bo  obtained  from  the  15  fringe  case  would 
be  to  the  nearest  fourth  of  a fringe  hut  since  noise  and  other  experimental 
errors  are  present,  only  about  1 fringe  resolution  can  be  obtained.  For 
the  case  of  5 fringes  in  the  projected  halo,  about  a quarter  of  a fringe 
resolution  is  possible.  The  system,  as  it  exists  in  its  present  state, 
can  analyze  data  at  a 2-3  point  per  second  rate  with  a very  low  resolution 
accuracy. 

Some  thought  has  gone  into  methods  of  improving  the  system  resolution 
which  appears  to  be  the  main  problem  at  the  present  time,  without  any  major 
change  in  the  overall  system  configuration. 

The  first  thought  which  came  to  mind  was  to  try  to  reduce  the 
secondary  speckle  noise  without  changing  the  fringe  resolution  of  the 
output  data.  Two  methods  of  in-house  experimentation  are  going  to  be 
performed  for  solving  this  problem.  The  first  will  be  to  slightly 
misdirect,  at  a high  fringe  rate,  the  readout  laser  beam  in  an  optical 
system  where  the  image  plane  is  held  stationary.  This  could  be  performed 
with  a motor  driven  rotating  wedge.  The  second  method  would  be  to  project 
the  output  fringe  pattern  onto  a liquid  crystal  dynamic  diffuser,  which 
acts  as  a random  phase  medium  that  is  electrically  controlled.  A surface 
of  this  type  would  tend  to  shift  the  secondary  speckle  a small  amount  in 
a random  fashion  at  a rapid  rate.  In  both  the  above  mentioned  methods, 
the  secondary  speckle  noise  would  be  minimized  by  taking  several  detector 
scan  data  samples  and  averaging  them  digitally  over  a short  period  of  time. 
This  process  could  be  done  with  a small  amount  of  modification  to  the 
existing  system. 

Another  way  of  analyzing  the  specklegram  data  that  will  be  considered, 
uses  the  optical  fourier  transform.  By  using  either  the  rotating  wedge 
or  the  diffuser  as  a noise  reduction  scheme,  an  optical  system  consisting 
of  a laser,  specklegram,  diffuser,  fourier  transform  lens,  D.C.  Block, 
and  appropriate  detector  could  perform  the  equivalency  of  the  digital 
fast  fourier  transform  operation,  in  an  optical  train  with  an  essentially 
zero  time  frame,  thus  eliminating  one  large  time  consuming  aspect  of  the 
present  method. 
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All  of  the  (ibove  techniques  are  beincj  considered  in-house  under 
project  2001-02-42  of  AFAL/OHO-2.  If  any  of  the  methods  prove  feasible  ' 

then  further  experimentation  could  be  performed  by  Flight  Dynamics 
Laboratory  approval. 
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A HASIRILITY  STUDY  OF  AUTOMATIC  FRINGF  COUNTING  TECHNIQUES 
1 3 November  1 973 
By  C.  Richard  Lane 


BACKGROUND: 

As  a result  of  structural  problems  associated  with  the  C-5  aircraft, 
Dr,  Lovelace  has  asked  that  all  laboratories  reflect  on  possible  solutions 
to  the  problem  of  monitoriny  the  structural  integrity  of  critical  areas 
on  the  aircraft.  Based  on  their  past  experience  in  the  area  of  non-destruc- 
tive testing  of  aircraft  structures.  Flight  Dynamic  Laboratory  organized 
an  approach  to  provide  such  a monitoring  system  in  cooperation  with  ARL, 
and  AFAL.  This  work  is  presently  being  done  under  AFFDL  work  unit  1467- 
02-12  with  Mr.  Gene  Maddux  and  Dr.  Frank  Adams  as  principal  investigators. 

The  initial  approach  taken  by  FDL  utilizes  the  recording  of  a double 
exposure  speckled  image  of  the  area  of  the  structure  being  monitored, 
where  some  strain  is  introduced  between  exposures.  This  double  exposed 
image  is  then  processed  to  make  a transparency.  Upon  illumination  by  a 
small  laser  beam  a set  of  fringes  are  projected  which  represent  the  corre- 
lation of  the  speckles  at  that  location  on  the  image. 

These  fringes  can  be  related  to  the  surface  motion  that  took  place 
between  the  two  exposures,  the  structural  integrity  at  that  point. 


Although  this  technique  is  limited  with  respect  to  the  size,  type, 
and  depth  of  the  flows  that  it  can  detect,  it  provides  a good  first  step 
and  is  essential  to  more  complex  techniques  that  have  potential  of  over- 
coming these  limitations. 
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For  each  C-5  aircraft  a large  number  of  these  specklegrams  would 
be  required  (5,000  - 10,000  if  one  specklegram  were  made  for  each  in- 
stalled wing  fastener)  and  a large  number  (100  - 1000)  of  points  per 
specklegram  would  need  to  be  analyzed.  It  is,  thus,  obvious  that  an 
automated  system  is  required  that  will  map  the  surface  motion  with  each 
specklegram  area. 

OBJECTIVE: 

The  objective  of  the  initial  in-house  effort  within  AFAL  was  to 
investigate  techniques  of  determining  the  spacing  and  direction  of  a 
projected  sample  fringe  pattern.  This  investigation  included  the  demon- 
stration of  the  techniques  on  a sample  specklegram  provided  by  AFFDL. 


APPROACH: 

As  a result  of  past  experience  with  similar  type  problems  in  the 
area  of  optical  processing  and  the  availability  of  equipment  within  AFAL 
two  approaches  to  the  problem  were  taken.  Both  approaches  required  that 
the  fringe  pattern  be  read-in  automatically  by  one  of  three  ways: 

a.  Scanning  point  detectors. 

b.  Electrically  scanned  linear  diode  arrays. 

c.  Electrically  scanned  2-dimensional  image,  dissector  tube. 

Once  the  fringe  information  was  detected,  two  approaches  were  con- 
sidered to  analyze  the  fringes  and  indicate  the  spacing  and  direction. 
The  two  approaches  were: 
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a.  Store  the  information  in  a memory  unit,  then  perform  a 2- 
dimensional  fast  fourier  transform  (FFT)  to  determine  the  spatial  fre- 
quency component  in  both  magnitude  and  angle.  The  FFT  was  chosen  over 
optical  transforming  due  to  the  lack  of  spatial  coherence  in  the  speckled 
fringe  pattern. 

b.  As  the  information  is  read-in,  perform  a peak  detection  opera- 
tion and  average  the  peak  spacing  for  scans  in  tv/o  directions.  By  means 
of  a simple  arithmetic  operation  the  maximum  spacing  and  direction  were 
then  determined. 

Evaluation  of  both  approaches  were  carried  out  on  the  sample  pro- 
vided. 


DISCUSSION  OF  THE  3 DEMONSTRATION  TECHNIQUES: 

Tv/o  of  the  three  demonstrations  to  be  discussed  now  were  done  only 
minimally  and  were  investigated  only  to  the  extent  of  being  considered 
as  a possible  solution  to  the  fringe  counting  problem.  Although  all 
three  demonstrations  were  capable  of  providing  a solution  to  the  problem, 
only  the  third  was  experimented  with  to  the  desired  extent.  It  will  also 
be  seen  that  the  3 demonstrations  were  similar  enough  in  theory  that  it 
can  be  assumed  that  if  one  experiment  works  then  all  three,  with  proper 
implementation,  can  serve  as  a solution.  All  three  techniques  have  been 
discussed  and  demonstrated  to:  Mr.  G.  Maddux  and  Dr.  Frank  Adams  of  the 
Flight  Dynamics  Laboratory. 

1)  On  6 August  1973  a sales  representative  from  Tektronix  Incorpora- 
ted accompanied  to  our  lab,  a processing  oscilloscope  being  corisidered 
for  purchase  under  project  2001-04-10.  The  sales  representative  set  up 
the  processing  system  in  the  lab  which  consisted  of  a Model  7704A 
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oscilloscope  with  standard  scope  plug-ins,  a Model  7001  Processor,  a POP 
’ 11-05,  16K  core  mini-computer,  a CRT  display  terminal,  and  a rapid  read- 

write  tape  punch.  The  self-scanning  detector  array  was  introduced  into 
the  system  through  the  oscilloscope  inputs.  The  array,  scanning  its  512 
elements  at  a 200  Hz  rate,  was  synchronized  to  the  scope  such  that  all 
■512  elements  were  displayed  on  the  scope  screen  in  one  sweep  (5m?ec). 

The  processor  portion  of  the  system  then  digitizes  the  data  displayed  on 
the  scope  screen  and  sends  it  to  the  POP  11  computer  along  with  a command 
of  which  software  program  to  implement  on  the  data.  The  resulting  ans- 
wer can  then  be  displayed  by  a curve  trace  on  the  scope  screen  or  CRT 
screen  of  the  terminal. 

Due  to  a lack  of  time  and  the  limited  software  programs  available 
by  the  sales  representative,  only  a fast  fourier  transform  program  was 
used  to  attack  the  fringe  counting  problem.  The  fringes  to  be  measured 
were  imaged  on  the  reticon  array  and  displayed  on  the  scope.  By  a stroke 
of  luck,  the  FFT  program  being  demonstrated  needed  512  digitized  points 
to  perform  the  transform,  exactly  the  number  of  detectors  in  the  scanning 
array.  The  signal  being  displayed  is  now  averaged  over  several  scans  of 
the  array  and  then  the  FFT  program  is  called  up  and  performed  on  the 
averaged  signal  and  the  result  displayed  with  512  points  on  the  oscillo- 
scope screen. 

Two  sets  of  fringes  were  experimented  upon  during  the  demonstration: 
An  equally  spaced  set  of  .05  inch  fringes;  and  a set  of  fringes  generated 
from  a speck  leg ram. 

a.  Figure  2a  shows  the  fringes  imaged  on  the  reticon  and  averaged 
1024  times  to  eliminate  unwanted  noise.  From  the  figure  there  are  21-22 
fringes  across  the  512  elements  of  data,  and  multiplying  the  number  of 
fringes  by  200  Hz  (the  scan  rate  of  the  detector),  gives  the  frequency 
that  the  computer  actually  sees. 
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(1)  21  X 200  = 4200  Hz  = 4.2  KHz 

(2)  22  X 200  = 4400  Hz  = 4.4  KHz 

Figure  2b  illustrates  the  fast  Fourier  transform  of  the  signal  of 
Figure  2a,  and  the  strong  peak  (arrow)  is  printed  out  as  4.3  KHz,  a 
reasonable  result. 

b.  Figure  3a  is  a scope  trace  of  a sample  specklegram  from  FDL, 
one  trace  a single  sweep  and  the  other  a trace  averaged  512  times.  These 
two  traces  illustrate  the  noise  reduction  capability  of  averaging  a sig- 
nal over  several  traces.  This  noise  reduction  factor  could  prove  very 
necessary  if  other  methods  of  fringe  counting  other  than  FFT,  is  imple- 
mented. The  FFT  displayed  in  Figure  3b  shows  two  significant  peaks  of 

data  relating  to  the  fringe  spacing.  The  peak  indicated  by  the  red  arrow 
is  a 200  Hz  frequency  spike  representative  of  the  scan  rate  of  the  reti- 
con  array.  This  frequency  spike  is  a consequence  of  the  input  fringes  not 
filling  the  entire  aperture  of  the  array,  thus  the  synchronizing  pulses 
from  the  clocking  circuit  of  the  array  are  prominent.  If  approximately 
the  first  30  reoisters  of  the  data  in  Figure  3b  are  set  to  zero  (200  Hz 
spike),  the  strongest  frequency  then  becomes  the  second  spike  and  this 
frequency  value  is  printed  in  Figure  3c  with  a 50mv  to  20  mv  vertical 
scale  change.  By  examination  of  Figure  3a  it  should  be  obvious  that  4.2 
MHz  represents  the  fringe  frequency  of  major  concern.  As  a simple  conclu- 
sion, a software  programming  addition  that  cancels  out  the  200  Hz  clock- 
ing rate  would  be  incorporated  and  only  the  desired  frequency  peak  would 
be  predominate.  By  retransforming  the  result  obtained  in  Figure  3c,  it 
is  seen  that  the  same  fringe  spacking  is  preserved  after  cancellation 
of  the  200  Hz  peak  (Figure  3d). 
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2)  The  second  method  of  demonstration  will  be  discussed  only  briefly 
since  it  is  similar  in  procedure  to  the  third  demonstration  with  slightly 
more  sophisticated  equipment.  Figure  lb  shows  the  basic  setup  used  for 
this  demonstration.  The  primary  difference  between  this  procedure  and  the 
first  demonstration  procedure  is  that  instead  of  calculating  a fast  fourier 
transform  with  a PDP  11  computer,  a counting  program  was  implemented. 

The  program  performs  an  adjacent  register  comparison  to  locate  peak  points 
and  then  checks  the  difference  between  these  peak  points.  This  approach 
is  instrumented  by  two  orthogonal  scans  across  the  fringe  pattern  with  an 
electronically  raster  scanned  image  dissector  tube.  Essentially,  the  re- 
sult is  an  X and  y fringe  spacing  measurement  of  which  an  actual  fringe 
spacing  may  be  calculated  along  with  a fringe  angle  displacement.  This 
procedure  will  be  discussed  more  thoroughly  in  the  next  demonstration. 

The  sole  purpose  of  this  experiment  was  to  indicate  the  feasibility  of 
using  an  image  disector  tube  for  the  fringe  counting  problem. 


3)  The  demonstration  to  be  discussed  now  utilized  equipment  already 
existing  in  the  laboratory  of  AFAL/TEL  and  thus  was  readily  available 
during  a longer  period  of  time  required  for  a more  thorough  investigation. 
The  equipment  used  and  its  configuration  are  illustrated  in  Figure  Ic. 

SYSTEM  DESCRIPTION: 

A laser  is  used  to  illuminate  a set  of  fringes  at  the  input  plane 
and  an  imaging  lens  projects  the  fringes  onto  an  x-y  scanning  aperture. 

The  variable  size  aperture  is  scanned  by  two  bi-directional  variable 
speed  D.C.  motors  which  are  n.anually  started  and  stopped.  The  aperture 
sizes  are  configured  such  that  the  resolution  of  the  system  can  vary  front 
a small  fraction  of  one  fringe  down  to  a single  fringe  or  essentially  no 
resolution.  This  aperture  size  variation  can  be  used  constructively  to 
reduce  noise  in  the  system  resulting  in  less  software  work  and  time  for 
the  computer  to  implement  at  a later  stage  in  the  processing.  A digital 


66 


AFAL-TR-77-203 


SCANNING  ELETRONICS 

I TELEVISION 
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Figure  lb 
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Y - SCAN  (MOTOR  DRIVE) 


PLOTTER  HP  9820 

CALCULATOR 


Figure  Ic 
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voltmeter  is  used  to  pick  up  the  intensity  pattern  of  the  fringes  and 
also  serves  as  an  analog-to-digital  converter  that  directly  interfaces 
with  the  Model  9820  Hewlett  Packard  Calculator.  The  digital  voltmeter 
and  calculator  combination  is  extremely  slow  for  use  as  a data  collector 
since  it  operates  at  a maximum  data  sampling  rate  of  3 samples  per  se- 
cond. This  fact  rules  out  this  type  of  system  as  a final  product  to 
solve  the  fringe  counting  problem,  although  it  strongly  demonstrates  the 
feasibility  of  using  a similar  fringe  counting  technique  with  a data 
reading  capability  of  approximately  4-5  magnitudes  greater.  A plotter 
can  also  be  tied  directly  to  the  system  giving  a visual  description  of 
the  fringe  data  that  the  calculator  will  process. 

SAMPLE  PROCEDURE: 

A good  method  of  examining  the  procedure  used  to  count  the  fringes 
would  be  to  look  at  Figure  4a,  b,  and  c,  which  list  the  calculator  program, 
in  some  detail. 

The  first  step  is  to  image  the  fringe  pattern  onto  the  scanning 
aperture  plane  and  position  the  aperture  in  the  y direction  ready  to 
scan  across  the  complete  fringe  pattern. 

By  simultaneously  running  the  program  and  starting  the  y direction 
scan  (manually),  the  "reads  data"  portion  of  the  program  begins  reading 
in  the  fringe  information  to  the  designated  registers  RO  - through  RIOO. 
Following  the  25  second  data  read  in  time,  the  program  initializes  step 

II  which  normalizes  the  data  such  that  the  minimum  voltage  register  is 
zero  and  all  other  registers  are  shifted  about  the  zero  level.  Section 

III  in  the  program  finds  the  register  with  maximum  value  and  affixes  this 
value  to  the  scale  setting  needed  for  plotting  the  100  data  registers  (IV) 
representing  the  fringe  pattern.  The  fifth  portion  points  out  the  data 
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0; 

0-X;FXD  0F 
1 : 

0-RX;  IF  (X+1>X)< 
249;  JMP  0F 
2:  ~ 

0^XF 
3: 

RED  3,A,RX;  IF  (X 
+ 1 >X)^100;  JMP  0F 
4:  n 

0 •X;RX  ►B;0-Ch 
5: 

IF  RX2B;RX-B;X^C 
F 

6: 

IF  (X+1 -X)<100; 

JMP  -Ir 
7: 

0-XF 

8; 

RX-B-RX;  IF  (X+1  ► 
X)'100;  JMP  0F 
9r  = 

0-X;  JMP  9f 
10: 

IF  RX>R130;RX 'R1 
30h 
11 : 

IF  {X+1-*X)<100; 

JMP  -IF 
12: 

SCL  0,  100,0,R130 
;0-XF 

13:  = 

PLT  X,RX;IF  (X+1 
-X)<100;  JMP  0F 
14:" 

PEN  ;SPC  2F 
15:  — 

0>X;FXD  0F 
16: 

IF  RX-.5E-2;PRT 
X;  JMP  IF 
17: 

IF  (X+1  •X)_100; 

JMP  -IF 

18:  

SPC  2;PRT  "THE  PEAKS 
ARE";PRT  "AT 
REGISTERS"F 

19:  

0-Xf 
20: 

IF  RX-  1E-2;0*RXf 
21 : " 

IF  lX+1  •X)'-100; 

JMP  -IF 


FIGURE  4a. 


•READS  DATA  FROM  DIGITAL  VOLTMETER 


II  VnORMALIZES  the  DATA  TO  ZERO 


III  >SETS  THE  SCALE  FOR  THE  PLOT 


IV  >PLOTS  THE  FRINGES 


^PRINTS  THE  WIDTH  OF  THE  FRINGE 
PATTERN  BY  REGISTERS 


VI 


^SETS  SMALL  VALUES  OF  NOISE 
AMPLITUDE  TO  ZERO 
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30: 

0-R(C+l);SPC  3; 
PRT  "DISTANCE  BE 
TWEEN" ;PRT  "ERIN 
GES  IN"h 


31 : 

ENT  "X-DIRECTION 
?",  B;  IF  B=l; 
PRT  "X-OIRECTION 
IS";JMP  2F 
32: 

PRT  "Y-DIRECTION 


IX  ASKS  FOR  DIRECTION  OF  SCAN 


CALCULATES  DIFFERENCE 
BETWEEN  PEAKS 


IF  SMALL  PEAKS  SHOW  UP, 
THESE  STEPS  SET 
SMALL  DIFFERENCES  TO 
ADJACENT  DIFFERENCE 


41 : 

200 -CF 
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► 


IF  B=1  ;JMP  2F 
45:  — 

0>R250;FXD  3; 

PRT  .018(R(A-2)/ 
(A-201 ))  ■R250; 

PRT"  INCHES"; 
SPC  8;  JMP  Z^■ 

46: 

0-R251;FXD  3; 

PRT  .029(R(A-2)/ 
(A-201 ))  ‘R251  ; 

PRT"  INCHES"; 

SPC  8F  _ 

47: 

STP  F 

48:  — 

PRT  "DISTANCE  BE 
TWEEN" ;PRT  "FRIN 
GES  IS"F 
49: 

PRT  R251/r(l+R25 

1R251/R250R250)* 

R252F 

50: 

PRT"  INCHES": 

SPC  I; PRT  "AN 
GLE  WITH";PRT  "R 
ESPECT  TO"; PRT" 

THE  X-AXIS"F 
51  : 

PRT  ACS  (R252/R2 
51);PRT" 

DEGREES" ;SPC  8F 
52;  — 


XIII 


XIV 


END  F 
R276 


FIGURE  4c. 

CALCULATES  SUM  OF  PEAK 
DIFFERENCES  AND  STORES 


CALCULATES  AND  PRINTS 
FRINGE  SPACING  IN 
BOTH  THE  X & Y DIRECTION 


from  X & Y fringe  SPACING 
THE  ACTUAL  FRINGE  SPACING 
AND  DIRECTION  ARE  PRINTED 
IN  INCHES  AND  DEGREES, 
RESPECTIVELY 
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registers  that  contain  fringe  data  or,  in  other  words,  registers  with  data 
-2 

more  than  5x10  in  amplitude.  This  information  merely  tells  the  program- 
mer if  noise  contained  in  the  systems  will  be  above  a certain  threshold  va- 
lue  set  by  the  prograniiier  (5x10  in  this  case),  and  immediately  following 
is  Section  VI  which  sets  these  noise  values  to  absolute  zero.  The  plots  and 
printouts  in  Sections  III,  IV,  AND  V are  mainly  used  to  check  the  program 
for  proper  operation  and  are  nonnally  eliminated  for  general  purposes  by  the 
jump  9 command  in  statement  9 of  Section  III.  The  first  printout  under  nor- 
mal operation,  performed  in  Section  VII,  calculates  the  registers  where  the 
peaks  of  the  fringes  occur  by  a method  of  adjacent  register  comparison,  and 
a resulting  printout  is  executed  in  Section  VIII.  An  example  of  this  print- 
out is  shown  in  Figures  5a  and  5b.  These  fringes  will  be  discussed  in  more 
detail  later  in  the  report.  Section  Ix  of  the  computer  program  simply  asks 
the  operator  the  direction  in  which  the  data  was  scanned  in  from  Section  1, 
and  correspondingly  enters  a scale  factor  that  has  been  manually  precalcu- 
lated to  convert  the  scanning  motor  speeds  into  a recognizable  format  such 
as  inches  for  this  case.  These  two  scale  factors  show  up  in  Section  XIII 
of  the  program  as  .018  for  the  y-direction  and  .029  for  the  x-direction.  In 
Section  X the  data  stored  in  Section  VII  are  substracted  and  these  differ- 
ences are  stored  again  for  later  use  in  Sections  XI  and  XII.  The  eleventh 
program  section  is  a safety  feature  installed  to  ensure  that  any  peaks  appear- 
ing as  a result  of  noise  in  the  system  that  were  not  eliminated  previously 
are  thrown  out  to  prevent  a false  representation  of  fringe  count.  The  fringe 
differences  are  now  averaged  over  the  entire  data  file  in  Section  XII  and  the 
uni-directional  (y)  spacing  is  calculated  with  the  motor  speed  scale  factor 
employed  to  give  an  answer  in  inches.  After  storing  this  spacing  in  inches 
the  entire  procedure  is  initiated  again  the  x-direction,  by  scanning  the  mo- 
tor drive  and  this  time  using  a different  scale  factor.  As  a result  of  Sec- 
tion XIII  both  X and  y fringe  spacing  information  is  stored  in  R251  and  R250, 
respectively.  Stepping  to  the  final  section  (XIV)  registers  R250  and  R251 
containing  the  x and  y information,  are  inserted  into  a simple  trigonometry 
function  to  determine  the  fringe  spacing  d and  fringe  angle  B which  are  the 
two  answers  being  sought  (see  Figure  4D). 
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Figure  4D 


Figures  5a  and  5b  illustrate  the  format  in  which  the  results  were  dis- 
played. The  answers  of  5a  were  made  from  fringes  spaced  at  .05  inch  apart 
and  5b  .10  inch  apart,  with  angles  for  the  two  sets  of  fringes  set  approxi- 
mately the  same. 

Finally  6 illustrates  a method  for  checking  the  fringe  angles,  by  sim- 
ply plotting  the  fringe  pattern  in  the  x and  y direction  and  graphically 
measuring  the  angle  offset. 
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CONCLUSION: 

Our  investigation  as  described  here,  indicates  that  several  techni- 
ques may  readily  be  applied  to  the  type  of  fringe  counting  problem 
considered. 

Several  choices  exist  for  reading  out  and  for  processing  ("counting") 
fringe  spacing.  The  choices  of  techniques  (and  optimum  hardware  implemen- 
tation of  same)  depend  heavily  on  the  detailed  technical  and  economic  para- 
meters and  constraints  of  the  application.  Consideration  of  these  details 
was  beyond  the  scope  of  these  initial  feasibility  investigations,  since 
all  equipment  used  in  the  demonstrations  was  presently  available  in-house. 

The  next  logical  step  appears  to  be  a design  analysis  to  determine 
the  optimum  system  configuration  taking  into  account  the  detailed  techni- 
cal and  economic  factors.  Following  this,  construction  and  testing  of  a 
prototype  system  would  be  appropriate. 

AFAL/TEL  would  be  pleased  to  undertake  further  work  on  this  problem 
along  the  lines  of  the  previous  paragraph  or  as  considered  appropriate  by 
mutual  agreement  of  AFAL  and  FDL. 
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APPENDIX  II 

CONTROL  PROGRAM  PRINTOUT 
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APPENDIX  III 

SYSTEM  OPERATING  INSTRUCTION 

OPERATING  PROCEDURES; 

I.  The  initial  set-up  should  be  performed  as  follows; 

1.  A 15  milliwatt  or  greater  Ue-Ne  laser  should  be  aligned  on  the 
centerline  of  an  optical  bench,  8.75  inch  above  the  surface  of  the 
bench,  with  all  of  the  optics  centered  on  this  beam. 

2.  Set  up  the  beam-reducing  optics  as  described  in  the  Optical 
Setup  Section  of  this  report  and  insert  an  optical  attenuator  imme- 
diately following  these  optics. 

3.  Place  the  x-y  stage  as  close  to  the  attenuator  as  needed. 

4.  Align  a ground  glass  screen,  with  a D.C.  spot  block,  so  that  the 
D.C.  component  of  the  fringe  halo  image  falls  on  the  block.  Position 
the  ground  glass  in  the  optical  axis  direction  so  that  a fringe  halo 
of  2-4  inches  in  diameter  exists. 

5.  An  imaging  lens  (f=28mm)  should  be  centered  on  the  optica?  axis 
and  positioned  to  image  the  fringe  pattern  onto  the  detectors.  The 
distances  between  the  lens  and  the  detector  all  depend  upon  the  image 
aperture  of  the  specklegram  recording  system.  Essentially  for  the 
fringe  spacing  of  a given  specklegram  to  be  analyzed  accurately,  the 
readout  halo  size  depends  on  the  recording  aperture. 

6.  Connect  the  x-y  stage  cables  to  the  appropriate  x-y  Aerotek  Drive 
Units  via  the  8 pin  connector  plugs. 

7.  A labeled  BNC  wiring  panel  is  attached  to  the  top  back  of  the  rack 
mounting  cart  and  is  illustrated  below. 


■0 

0 

<o> 

0. 

0 

x-direction  x-clock 

y-direction 

y-clock 

ret  Icon  clock 

0 

<o) 

0 

x-video 

y-vldeo 

channel  2 
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A list  of  to  and  from  BMC  connections  are  shown  in  the  followiiHj  table; 


WIRING  DIAGRAM 

FROM 

TO 

Aerotek  Drive-(x)  DIR 

Aerotek  Drive-(x)  CL 

Aerotek  Drive- (y)  DIR 

Aerotek  Drive-(y)  CL 

Reticon  -x  video 

Reticon  -y  video 

Reticon  -Ext.  clock 

Reticon  -EOS 

Panel  - x - direction 

Panel  - x - clock 

Panel  - y - direction 

Panel  - y - clock 

Panel  - x - video 

Panel  - y - video 

Panel  - reticon  clock 

Panel  - channel  2 

B.  AC  connections  from  the  Aerotek  Drive  Units,  PDP  11  computer, 
off-on  indicator  light,  and  control  electronics  panels  should  be 
plugged  into  the  AC  terminal  strip  in  the  lower  rear  of  the  rack 
mounting  cart. 

9.  Connect  the  teletype  I/O  to  the  computer  I/O  cable  and  plug  the 
teletype  into  a wall  socket. 

10.  Turn  on  all  equipment  and  check  out  this  equipment  with  a test 
paper  tape  given  to  you  by  AFAL/DHO-2. 


II.  AFTER  CHECK-OUT  (Initializing  for  a read  sequence) 

1.  Place  Specklegram  into  4 inch  by  5 inch  holder. 

2.  Align  the  readout  beam  ti,  the  upper  right  hand  corner  of  the 
specklegram  plate,  (facing  the  plate),  by  setting  the  local-remote 
switch  on  the  x drive  unit  to  the  local  position  and  pushing  the  slew 
button  until  the  upper  right  hand  x position  is  reached.  The  + or  - 
switch  will  be  set  according  to  the  initial  position  of  the  plate 
before  setting  up.  The  slew  speed  can  also  be  changed  by  the  small 
potentiometer  on  the  front  panel.  Repeat  for  the  y-position. 

3.  Reset  the  x and  y digital  display  to  zero  by  depressing  the  reset 
button  and  switch  the  remote-local  switches  to  the  remote  position. 
The  system  is  now  ready  for  computer  control  operation. 
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4.  Load  the  control  program  through  the  paper  tape  reader  on  the 
teletype  by  the  following  procedure: 

A. )  If  the  Bootstrap  loader  Is  not  contained  in  core-memory, 
a hand  loading  procedure  is  as  follows; 

a. )  Toggle  in  017744  onto  the  computer  switches  and  hit 
LOAD  ADDR  Button. 

b. )  Toggle  016701  and  hit  the  DEP  button  and  continue  DEP 
operations  on  these  numbers:  000026,  012702,  000352,  005211, 
105711,  100376,  116162,  000002,  017400,  005267,  177756,  000765, 
177560,  000000. 

B. )  Set  the  ABSOLUTE  LOADER  paper  tape  into  the  reader  with  the 
leader  section  in  the  reader  section,  and  set  the  reader  switch 
to  the  top  position  called  START. 

C.  ) 017744  and  LOAD  ADDR  then  START.  The  Absolute  loader  should 

now  pass  through  the  reader  and  stop  when  finished. 

D. )  Load  the  Control  tape  into  the  reader  with  the  leader  posi- 
tioned in  the  reader  section,  and  set  the  read  switch  up. 

E.  ) 017500  and  LOAD  ADDR  then  START.  The  control  tape  should  now 

pass  through  the  reader  and  stop  when  finished.  The  control  program 
is  now  loaded. 

5.  Make  sure  all  connections  are  made  and  all  power  switches  are 
activated.  The  program  is  now  ready  to  run. 

6.  000500  and  LOAD  ADDR,  then  START. 

7.  The  teletype  will  now  ask  initializing  questions  before  displaying 
specklegram  data  and  these  questions  can  be  answered  as  described  in  the 
Software  Section  of  this  report.  When  the  data  run  is  finished,  the 
teletype  will  print  "SCAN  COMPLETED". 

8.  If,  before  shutdown,  another  plate  is  to  be  analyzed,  only  the  fol- 
lowing steps  have  to  be  repeated; 

II.  - 1,2, 3. 6 

9.  The  computer  can  be  stopped  during  any  phase  of  operation  by  simply 
hitting  the  HALT  Button  and  after  a data  run  there  is  no  specific  shutdown 
sequence  except  that  the  HALT  button  should  be  activated  before  shutting 
any  other  power  off. 
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10.  Once  the  procedures  of  Sections  I and  II  of  the  operating  pro 
cedures  have  been  performed  complete,  only  steps  II.  - 1,  2,  3,  4B,  4C, 
4D,  4E,  5,  and  6 have  to  be  repeated. 
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APPENDIX  IV 

A TYPICAL  SPECKLEGRAM  EXPOSURE  AND  READOUT 

Appendix  V is  designed  to  inform  the  reader  in  brief  of  a typical 
experiment  run  through  from  initial  exposure  to  final  readout,  w th  the 
understanding  that  specific  details  can  be  obtained  from  the  references 
listed  at  the  end  of  the  report. 

1.  EXPOSURE 


In  recording  speckelgram  data,  if  the  laser  has  a wavelength  of 
5140  angstroms  and  the  lens  has  an  f-number  equal  to  2,  then  the  speckle 
size  projected  onto  the  film  plane  is  approximately  1 micrometer,  which 
is  defined  as  the  minimum  speckle  size  for  our  purposes.  This  minimum 
sets  one  paraneter  related  to  the  system  resolution.  The  next  question 
related  to  the  system  resolution  is  what  is  the  maximum  and  minimum  shift 
that  can  be  tolerated?  Sampling  Theory  requires  that  the  shift  be  greater 
than  2 speckle  sizes  (2  micrometers).  The  maximum  shift  is  defined  by  re- 
quiring the  shift  not  to  be  larger  than  the  readout  beam  size  (.5  to  1 
millimeter).  The  reason  for  this  maximum  is  that  the  surface  across  the 
reading  aperture  can  be  shifted  in  only  one  direction  with  the  randomness 
of  the  speckle  remaining  unchanged.  The  exposure  time  required  on  the 
film  plane  is  proportional  to  the  film  sensitivity  and  for  a typical  film 
(Kodak  649-F),  500  ergs/cm  is  needed.  Assuming  now  for  simplicity  a 
square  centimeter  of  gray  steel  (55*  reflectivity)  as  a specimen,  a 50 

milliwatt  laser  as  an  illuminator  source,  in  which  all  of  the  light  inter- 
2 

cepts  the  1 cm  area,  and  a cosine  reflection  distribution  from  the 
specimen  with  1*  of  the  reflected  light  reaching  the  imaging  lens,  a 
typical  exposure  might  be  as  follows: 

-4 

1000  ergs  = 10  watt-sec 
500  ergs  = .05  milliwatt-sec 
.05/50/. 55/. 01  = .18  seconds 

no 


i 

i 
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The  .18  second  exposure  time  assumes  ideal  conditions,  thus  optical 
losses  in  the  system  would  require  more  exposure  time. 


The  processing  required  is  the  standard  film  process: 

1.  Develop  (D-19)  - 6 minutes 

2.  Stop  (Stop  Bath)  -30  seconds 

3.  Fix  (Rapid  Fixer)  - 5 minutes 

The  plate  is  now  ready  for  readout. 

2.  READOUT 

Readout  techniques  used  for  extracting  the  fringe  information  from 
a speckl egram  and  a typical  fringe  halo  are,  for  the  purposes  of  this 
report,  as  follows.  A 1 mm  readout  beam  diameter  is  chosen,  thus  providing 
a maximum  of  13000  data  points  on  one  4 inch  by  5 inch  film  plate.  Since 
the  readout  procedure  for  each  point  is  the  same,  the  procedure  for  only 
one  point  will  be  discussed.  The  halo  containing  the  fringe  information 
is  defined  as  a conical  ray  of  light  diverging  from  the  center  of  the  read- 
out beam  with  a cone  angle  directly  calculatable  from  the  recorded  speckle 
size  as  follows: 

2>/D  = 0 

where  D=the  speckle  size.  If  D=1  micrometer  and  X = .5  micrometer, 
then  >1  = 1 radian. 

Since  the  recording  speckles  are  arranged  in  a random  fashion  and 
the  speckle  shift  is  in  only  one  direction  over  the  area  defined  by  a 1 
mm  aperture,  then  the  analogy  can  be  made  that  the  random  speckle  shift 
acts  similarly  to  a double  slit  optical  configuration. 


The  equation  for  the  double  slit  is 

n ' = 2d  s i n 6 
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where  n is  the  order  of  the  fringe,  ).  is  the  wavelength  of  light,  d is 
the  distance  between  slits  (or  speckle  correlation  shift)  and  9 the  angle 
of  the  fringe  transmission. 

This  says  that  the  first  dark  fringe  should  be  located  at  an  angle 
of 


9 = A/2d 


if  9 is  a small  angle. 


It  should  be  obvious  at  this  point  that  the  specklegram  strain  analy- 
sis technique  is  an  excellent  tool  for  looking  at  in-plarie  displacements 
or  strains  parrallel  to  the  recording  surface,  while  holographic  interfer- 
ometry best  applied  to  out-of-plane  measurements;  thus  the  two  techniques 
should  complement  each  other  for  solving  many  strain  and  stress  problems 
existent  in  the  Air  Force. 
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